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ABSTRACT 

in 

176P/LINEAR is a member of the new cometary class known as main-belt comets (MBCs). It displayed cometary activity shortly 
during its 2005 perihelion passage, which may be driven by the sublimation of subsurface ices. We have therefore searched for 
emission of the H 2 1 io — lot ground state rotational line at 557 GHz toward 176P/LINEAR with the Heterodyne Instrument for the 
Far Infrared (HIFI) onboard the Herschel Space Observatory on UT 8.78 August 201 1, about 40 days after its most recent perihelion 
passage, when the object was at a heliocentric distance of 2.58 AU. No H,0 line emission was detected in our observations, from which 
we derive sensitive 3-cr upper limits for the water production rate and column density of < 4 x 10 25 mol s -1 and of < 3 x 10'° cirr 2 , 
respectively. From the peak brightness measured during the object's active period in 2005, this upper limit is lower than predicted 
by the relation between production rates and visual magnitudes observed for a sample of comets at this heliocentric distance. Thus, 
176P/LINEAR was most likely less active at the time of our observation than during its previous perihelion passage. The retrieved 
upper limit is lower than most values derived for the H 2 production rate from the spectroscopic search for CN emission in MBCs. 
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1. Introduction Theoretical models suggest that the snow line was initially 

close to the Mars orbit due to the absorption of stellar radiation 

^ Classical main-belt asteroids are small bodies that orbit the Sun by dust Kasselov & Lecar|2000l iLecar et al.|2006b. Therefore, 

^ in low inclination and low eccentricity orbits between the or- objects formed at their current locations in the outer asteroid belt 

\l bits of Mars and Jupiter. Physically, asteroids are thought to be may haye been able tQ accumulate water ice in subs urface reser- 

. devoid of volatiles, while comets are icy bodies that become yoirSi despite the effect of solar radiation . Determining the com- 

OO active in the inner solar system thanks to the sublimation of position of tMs clasfj of objects can proyide important clues to 

^ ices, mostly water. Comets originate in the outskirts of the so- both the thermal properties that allow water t0 survive in subsur . 

1_I lar s y stem be y° nd the snow llne ' where temperatures in the solar face j and the distl - ibution of volatile materials in the solar 

. . nebula were low enough for water to condense onto icy grains nebula tQ constrain planet formation mechanisms. Additionally, 

> ( |Hayashi|1981^ A new class of bodies has been discovered re- MBCfj may haye played an important role in the delivery of wa . 

ITZ cently, the so-called main-belt comets (MBCs), which have or- ter and other volatiles t0 the inner solar systenii inc i ud i ng the 

rN bital properties that are indistinguishable from standard asteroids E^th 
~3 with a Tisserand parameter with respect to Jupiter that is greater 

than three, and they display cometary activity in the form of a The MBC 176P/LINEAR (hereafter 176P) was discovered 

dust tail during pail of their orbit. Numerical simulations have ln 1999 and originally categorized as asteroid 1 18401 LINEAR, 

shown that these objects are not comets from the Kuiper Belt or This ob J ect belongs to the Themis asteroid family. Cometaiy 

Oort Cloud that have been recently transferred to orbits within activity was report ed for this object around perihelio n in 2005 

the main belt, but instead are most likely formed in situ at their (Hsi eh et al.|20Tl) by the Hawaii Trails project (HTP; |Hsieh & 



current locations ( [Fernandez et al.|2002] >. |Jewitt||20 06 1 Hsie h||2009| ). It displayed a mean photometric ex 

cess of ~ 30% during a month-long active period around its per- 
* Herschel is an ESA space observatory with science instruments ihelion passage, consistent with an approximate total dust mass- 
provided by European-led Principal Investigator consortia and with im- loss of ~ 7 X 10 4 kg ( |Hsieh et al.|201 1| >. Although ice sublima- 
portant participation from NASA. tion is expected to trigger MBC activity, gas emission has never 
** Current address: Department of Astrophysical Sciences, Princeton been directly detected in these objects owing to their low activ- 
University, NJ 08544, USA ity, which requires very sensitive observations. Herschel proves 
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to be the most sensitive instrument for directly observing water 
in a distant comet (e.g. Bockelee -Morvan et al.|[20T0| . In this 
paper we present the Herschel observation of the 1 10— lot fun- 
damental rotational transition of H 2 at 557 GHz in 176P. This 
observation is intended to test the prediction that the observed 
cometary activity of MBCs is driven by sublimation of water 
ices and to constrain the production process. 

2. Observations 

The MBC 176P was observed with the Heterodyne Instrument 
for the Far Infrared (HIFI; |de Graauw eTaLll^OlOl ), one of the 
three instruments onboard the ESA Herschel Space Observatory 
(Pilbr att et al.||20T0] l, within the framework of the Herschel 
guaranteed-time key program "Water and related chemistry in 
the solar system" (Hartogh et al. 2009]). HIFI provides very 
high-resolution spectroscopy that can resolve the line shape 
and enable the determination of accurate production rates (e.g., 
|Hartogh et al.|2010| l. 176P was the best MBC target in terms of 
its visibility close to the perihelion passage and anticipated line 
strength to be observed by Herschel during the mission lifetime. 
It passed its perihelion on 30 June 201 1 at a distance of 2.57 AU 
from the Sun and was observed by Herschel on UT 8.78 August 
201 1 with a total on-target integration time of 4.8 hours, when it 
was at a heliocentric distance of 2.58 AU and a distance of 2.55 
AU from the satellite (Herschel ObsID 1342225905). The object 
was tracked using an up-to-date ephemeris provided by the JPL's 
Horizons system. 

The line emission from the fundamental (ljo-loi) rotational 
transition of ortho-water at 557 GHz was searched for in the up- 
per sideband of the HIFI band la mixer. The observation was 
performed in the frequency-switching observing mode with a 
frequency throw of 94.5 MHz, using both the wide band spec- 
trometer (WBS) and the high-resolution spectrometer (HRS). In 
this observing mode there is no need to observe a reference posi- 
tion on the sky and the on-target integration time is maximized. 
However, the statistical noise may be underestimated for ob- 
servations in frequency-switched mode owing to uncertainties 
in baseline removal (Bockelee-Morva n et al.||2012| l. The spec- 
tral resolution of the WBS is 1 MHz (~ 0.54 kms _1 at the fre- 
quency of the observed line), while the HRS was used in its 
high-resolution mode with a resolution of 120 kHz (~ 0.065 
kms -1 ). The main beam brightness temperature scale was com- 
puted using a beam efficiency of 0.75 and a forward efficiency 
of 0.96. The folded spectrum was obtained by averaging the 
original spectrum with a shifted and inverted copy. Horizontal 
and vertical polarizations were averaged, weighted by the root 
mean square amplitude, to increase the signal-to-noise ratio. The 
pointing offset of horizontal and vertical polarization spectra is 
6'.' 6 in band la, corresponding to approximately 20% of the half- 
power beam width at the observed frequency. 

3. Data analysis 

The data analysis was performed using the Herschel Interactive 
Processing Environment (HIPE) software package (Ott 20 1 0) > . 
The standard HIFI processing pipeline v7.1.0 was used to re- 
duce the data to calibrated level-2 data products. The frequency- 
switching observing mode introduces a strong baseline ripple. 
To obtain a reliable estimate of the noise present in the mea- 
sured data, the baseline has to be removed, which is usually ac- 
complished by fitting a linear combination of sine waves using 
the Lomb-Scargle periodogram technique. Nevertheless, the in- 
strumental processes responsible for the baseline are in general 
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Fig. 1. Original HRS subband 1 spectrum of the H 2 1 to— lot 
line at 556.936 GHz observed on UT 8.78 August (upper panel), 
and several low and high frequency components of the spectrum 
determined using the EMD analysis (four lower panels) with la- 
bels indicating the mode number. The vertical axis is the cali- 
brated main beam brightness temperature. The lower horizontal 
axis is the upper sideband frequency, while the upper axis shows 
the velocity with respect to the nucleus's rest frame. 



combinations of linear distortions of different components in the 
receiver/spectro meter subsystems with a small fraction of non- 
linear processes, which may cause an aperiodicity in the ripple. 
Analyzing such contaminated signals by assuming a linear re- 
lation among the signal components (a fundamental assumption 
for all the Fourier-based techniques) is not always suitable, de- 
pending on the degree of nonlinearity. In this work, we utilize a 
relatively novel approach specifically developed for analysis of 
aperiodic and nonlinear signals - the Hilbert-Huang Transform 
(HHT ; |Huang et al.|[T998l |Huang et al.|[T999t |Battista et all, 
2007). This approach combines the empirical mode decompo- 
sition (EMD) procedure, which decomposes the original signal 
into its intrinsic mode functions (IMFs; representing the differ- 
ent modes of oscillations) with the Hilbert transform that can be 
then used in computing the instantaneous frequencies. The EMD 
technique extracts all the oscillatory modes, including all the 
baseline ripple components (the smooth low-frequency modes), 
as well as the highest frequency components, usually the noise. 
Another important property of IMFs is that they obey a sim- 
ple additive rule to reconstruct the original signal exactly. This 
makes the EMD approach an accurate and versatile tool that 
has been successfully used in many areas from atmospheric sci- 
ence to cosmology such as denoising and detrending, as well as 
the time series analysis tool for identifying periodic and quasi- 
periodic features (see Duffy 2004 ; Batti sta et al.|2 007 , and refer- 
ences therein). 

In this work we applied the EMD technique to the measured 
WBS and HRS spectra to obtain the highest frequency IMFs 
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Fig. 2. High-frequency component of the HRS spectrum of the 
H 2 lio-loi line at 556.936 GHz observed on UT 8.78 August 
with overplotted synthetic spectrum of the 3-cr upper limit 
shown as the dashed line. The vertical axis is the calibrated main 
beam brightness temperature. The lower horizontal axis is the 
upper sideband frequency, while the upper axis shows the veloc- 
ity with respect to the nucleus's rest frame. 

(dominated by Gaussian noise), which are usually the first modes 
( jDuffy |2004| l . Some of the modes obtained from the decomposi- 
tion of the HRS spectrum are shown in the lower panels of Fig. [T] 
Then, the root mean square (rms) noise of the brightness temper- 
ature in each spectrum is used to derive a 3-cr upper limit for the 
H 2 production rate shown in Table [T] The brightness tempera- 
ture rms of the WBS and HRS spectraaiffer owing to their differ- 
ent spectral resolutions as predicted by the Herschel Observation 
Planning Tool (Hspot). The rms agrees with the value derived 
using the Lomb-Scargle periodogram method to determine the 
frequencies of the baseline. Typically the Lomb-Scargle imple- 
mentation requires 20 to 30 sinusoidal components to achieve a 
good fit of the baseline ripple, while the HHT analysis provides 
a good estimate of the noise with eight modes and has an addi- 
tional advantage in computing speed. Fitting a narrow frequency 
range of ~ 40 MHz around the water line with a high-order poly- 
nomial tends to underestimate the noise level determined from 
the whole band by about a factor of two. However, the differ- 
ence introduced by the baseline fitting method is smaller than 
the uncertainty in the upper limit for the outgassing rate due 
to unknown model parameters, and therefore our data reduction 
method does not modify the conclusions. The EMD and HHT re- 
duction methods applied to the baseline removal and denoising 
of the HerschelfHIFl data will be described in detail in a future 
work (Rezac et al. in preparation). We show the baseline sub- 
tracted HRS spectrum in Fig. [2] with the expected line emission 
overplotted. 

4. Outgassing rate 

There is no evidence of H z O emission in our observation, al- 
though it is expected that the object's dust emission activity is 
driven by the sublimation of subsurface material as it approaches 
perihelion. A molecular excitation model based on the publicly 



available accelerated Monte Carlo radiative transfer code ratran 
( jHogerheijde & van der Tak|2000| l is used to calculate the pop- 
ulation of the rotational levels of water as a function of the nu- 
cleocentric distance. The code includes collisional effects and 
infrared fluorescence by solar radiation to derive the produc- 
tion rates. We used the one-dimensional spherically symmet- 
ric version of the code following the description outlined in 
Bensch & Bergin ([2004) that has been used to analyze Herschel 
and ground-based cometary observations (see e.g. Hartogh et al. 
20101 |de Val-Borro et al.|[20T0l |Hartogh et aLl|2011[ |de VaT 
Borro et al.||2012| >. The model input parameters are the gas ki- 
netic temperature, which controls the molecular excitation in the 
collisional region, and the electron density. We assume a gas ki- 
netic temperature in the range 20^40 K. The electron density and 
temperature profiles from |Biver ( 1997 1 are adopted. Since the 
electron density in the coma is not well constrained, an electron 
density scaling factor of x„ e = 0.2 with respect to the standard 
profile derived from observations of comet lP/Halley has been 
used (e.g. |Hartogh e t al. 2010). The expansion velocity is as- 
sumed to be constant in the coma, and the radial gas density pro- 
file for water was obtained using the standard spherically sym- 
metric Haser distribution. The gas expansion velocity derived 
from the pre- and post-perihelion evolution of comet C/1995 Ol 
(Hale-Bopp) is given by v eX p = 1.12 x r h 0Al kms _1 (Biver et al. 
1997f l, which corresponds to 0.76 kms~'at 176P's heliocentric 
distance. Since Hale-Bopp was a very active comet, this expan- 
sion velocity is a very conservative upper limit. A thermal veloc- 
ity of 0.35 kms -1 is obtained from the temperature expected at 
the subsolar point where ice sublimates. Then, as cometary at- 
mospheres are formed by quasi-adiabatic expansion, the gas ac- 
celerates as it expands, as observed with the large field of view of 
the 18-cm OH observations. For low-activity and distant comets 
an expansion velocity close to 0.5 km s _1 is determined from the 
shapes of the OH l ine observed with the Nancay radio telescope 



( Tseng et al.|2007 1, but observations for Qqh < 10 28 molec. s 



are lacking, and observations at rj, > 2 AU are rare. Odin obser- 
vations of the H 2 557 GHz line toward the active comet C/2003 
K4 (LINEAR) at 2.2 AU from the Sun are co nsistent with an ex - 
pansion velocity on the order of 0.5 kms'l fBiver et aL|2007 i. 
The dust particles ejected by 176P have an approximate velocity 
or ~ 5 m s^calculated from numerical simulations to match op- 
tical observations (Hsie h et al.|201 l| l. However, dust is expected 
to be much slower than gas if it is formed by large and heavy 
particles with very low outgassing rates. 

The mean of the derived upper limits on the total H 2 col- 
umn density integrated within the beam for the considered model 
parameters is (Nu^o) < 3xl0 10 ctrT 2 . For a low expansion veloc- 
ity characteristic of weak comets, we derive a sensitive 3-cr up- 
per limit on the water production rate of < 2.1 x 10 25 molec. s _1 
from the WBS and HRS data (see Table [T}. An upper limit of 
< 4 x 10 25 molec. s _1 is derived from the mean of the WBS and 
HRS upper limits for gas expansion velocities between 0.4-0.7 
km s 1 and gas kinetic temperatures between 20^-0 K. 

With the exception of 133P/Elst-Pizarro, our upper limit for 
the H 2 production rate is more stringent than any of those 
derived in other MBCs from the spectroscopic search for CN 
emission in the optical when the objects were active. In addi- 
tion, it does not require the uncertain assumption of a 2cn/2h,o 
value. Adopting a 2cn/6h,o mixing ratio of 10~ 3 the upper 
limits are g H ,o < 1-3 X 10 24 molec. s" 1 in 133P/Elst-Pizarro 
( jLicandro et al.|201l}, Q H ,o < 1-4 x 10 26 molec. s 1 in P/2008 



Rl (Garradd)(Jewittetal.2009i, <2h o < 1.3xl0 27 molec. s _1 in 
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Table 1. Standard deviation of the brightness temperature and 
line area, and retrieved 3-<x upper limits of the H 2 (D production 
rate in 176P/LINEAR. 



Spectrometer 




°7 T mB dv 


2h 2 o 






(K) 


(Kkms" 1 ) 


(mol s" 


- 1 ) 


WBS 


6.598 x lfr 4 


6.172 x 1(T 4 


< 2.08 x 


10 25 


HRS 


1.998 x 1(T 3 


6.365 x 10~ 4 


< 2.14 X 


10 25 



Notes. (fl) Production rates derived for a gas kinetic temperature of 20 
K, expansion velocity of 0.5 kms~', and an electron density scaling 
factor of x„, = 0.2. 



2012b), £>h,o < 6x 10 26 molec. s 1 
^ 2012c), and Q H ,o < 



P/2006 VW139 (jHsieh et al._ 
in P/2010 R2 (La Sagra) ( jHsieh et al 

9 x 10 26 molec. s" 1 in the collisionally disrupted main belt object 
(596) Scheila (Hsieh et a l.|2012a) . However, no direct search for 
H 2 has been carried out before, and there are substantial un- 
certainties in the estimation of the water production from the 
CN emission given the wide range of observed 2cn/Gh,o ra- 
tios in comets and their dependence on heliocentric distance. 
Therefore, these previously published values can only be con- 
sidered as order-of-magnitude approximations. 



5. Discussion 

We observed water emission in 176P with Herschel/HWl to test 
the prediction that cometary activity in MBCs is driven by sub- 
limation of water ices from the nucleus. There are several mech- 
anisms that have been proposed to drive mass loss from small 
bodies, including sublimation of subsurface ices, rotational in- 
stability, impact ejection and thermal fracture (see |Bertini|2011| 
|Jewitt|2072] for recent reviews of MBC physical properties and 
activation mechanisms). The cometary activity observed in 176P 
was initially found to suggest the presence of sublimating sub- 
surface ice that may have been exposed by recent collisions 
( Hsieh et al.||201 1 1. This view is supported by the detection of 
water ice absorption in spectroscopic observations centered at 
3.1 fim of the surface of the largest asteroid of the Themis as- 
teroid family, 24 Themis, which belongs to the same dynamical 
family as 176P ( |Rivkin & Emery|[20T51 |Campins et "^[2010) , 
although it has been claimed that the measured spectra are con- 



sistent with the transmission spectra of goethite ( |Beck et al. 



2011 



From the search for the H 2 Iiq- 



loi rotational line at 



557 GHz in 176P, a 3-cr upper limit for the H z O production rate 
of < 4 X 10 25 molec. s _1 is derived from the WBS and HRS data, 
for gas expansion velocities between 0.4—0.7 kms'and gas ki- 
netic temperatures between 20 and 40 K. Using the peak value 
of the /?-band magnitude, m(l, 1,0) = 15.35 + 0.05, measured 
when the object was active in late 2005 at a heliocentric dis- 
tance of 2.58 AU, a V-band magnitude of my(l,rh,0) = 17.8 
is obtained ( jHsieh et al.||2009) . Since the cometary activity in 
176P is indicative of ice sublimation with a 30% contribution 
of the coma to the total brightness, the scaling relation between 
gas production rates and heliocentric magnitudes from |Jorda| 
|et al.| ( |2008) predicts a water production rate of approximately 
1.0 x 10 /o molec. s _1 . This correlation has been obtained for a 
sample of 37 comets with heliocentric distances between 0.32- 
4.53 AU. Thus, if ice sublimation is the driving mechanism of 
176P's activity, the derived H 2 production rate is too low by 
about a factor of two to explain the activity level during its 
2005 perihelion passage. It is unlikely that sublimation of car- 
bon monoxide ice is the source of the activity in MBCs since 



the temperature in the region where the comet formed could 
not have allowed condensation of CO. We note that the IJordal 
|et al.| ( |2008] ) correlation should be taken with some care because 
none of the comet measurements were obtained at the bright- 
ness level of 17 6P. From the dus t production rate of 0.07 kg 
s _1 estimated by Hsieh et al. ( 201 1) , a water production rate of 
2.3 x 10 24 molec. s" 1 is inferred assuming a dust-to-gas ratio 
of one. However, there are significant uncertainties in the value 
of the dust production rate determined from photometric data, 
and additionally the dependence on the dust-to-gas ratio with 
heliocentric distance is poorly constrained. Mid-infrared photo- 
metric observations of 176P on 23-24 April 2010 by the Wide- 
field Infrared Survey Explorer (WISE) mission did not provide 
any indication of a coma ( Bauer et al.||2012 i, and there are no 
other published infrared or optical data closer to the last peri- 
helion passage. We conclude that water was not detected in our 
observation because the water production rate was lower than 
4 x 10 25 molec. s _1 or the object was not active during our ob- 
servation, so a more detailed study is needed to shed light on the 
activation mechanism in MBCs. 
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